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This is a follow-up to our earlier work (fil 0, for the energies and the charge (vector) and matter 
(scalar) distributions for S-wave states in a heavy-light meson, where the heavy quark is static and 
the light quark has a mass about that of the strange quark. We now study the radial distributions of 
higher angular momentum states, namely P- and D-wave states. In nature the closest equivalent 
of this heavy-light system is the B s meson. 

The calculation is carried out with dynamical fermions on a 16 3 x 32 lattice with a lattice spacing 
of about 0.10 fm generated with the non-perturbatively improved clover action. It is shown that 
several features of the energies and radial distributions are in qualitative agreement with what one 
expects from a simple one-body Dirac equation interpretation. 
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1. Energies 

The basic quantity for evaluating the energies of heavy-light mesons is the 2-point correlation 
function - see Fig. [T[ It is denned as 

C 2 (T) = (P ! rG q (x,t + T,t)P t+T r*U Q (x,t,t + T)), (1.1) 

where U@(x,t,t + T) is the heavy (infinite mass)-quark propagator and G q (x, t + T,t) the light anti- 
quark propagator. P t is a linear combination of products of gauge links at time t along paths P and 
r defines the spin structure of the operator. The (...) means the average over the whole lattice. The 
energies are then extracted by fitting the C2 with a sum of exponentials, 

w max 

c 2 (r)f* £ c ie - m < T Ci , (i.2) 
;=i 

where N max = 2-4, T < 14. 

Calculations were made using three different lattices with /3 = 5.2, C§y^ = 2.0171 non- 
perturbatively improved clover fermions. The parameters are given in Table [I] and the extracted 
energies are summarized in Fig. |2[ The notation L+(— ) means that the light quark spin couples 
to orbital angular momentum L giving the total j = Ldb 1/2. 2S is the first radially excited L= 
state. Energies are given with respect to the S-wave ground state. The plot also shows a comparison 
between the static (infinitely heavy) and smeared ("sum6") heavy quark. The "sum6" is APE type 
smearing, where the simple gauge links in the time direction are replaced by the sum of six gauge 
link staples. The maximum distance from the original, simple link is one lattice spacing. The effect 
of the smearing on energy differences seems to be small. 



(x,t + T) (X,i!) 




Figure 1: The two-point correlation function C2 (on the left) and the three-point correlation function C3 (on 
the right). The probe ® is at a distance R from the heavy quark. 
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DF3 


0.1350 


0.110 


1.1 


1.93(3) 


DF4 


0.1355 


0.104 


0.6 


1.48(3) 


DF5 


0.1358 


0.099 


0.3 


1.06(3) 



Table 1: Lattice parameters. We use ro = 0.525(25) fm to convert lattice results to physical units. 
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Figure 2: The energies of different angular momentum states with respect to the S-wave ground state. The 
error bars show the statistical error and the estimated systematic error from the extraction procedure. No 
continuum limit is taken. See Section |3| for more details about the Dirac model fit. 

2. Radial distributions 

For evaluating the radial distributions of the light quark a 3-point correlation function is needed 
- see Fig. [jj It is defined as 

c 3 (R,T) = (r^u Q rG ql eG q2 ). (2.1) 

We now have two light quark propagators, G q \ and G q 2, and a probe &(R) at distance R from the 
static quark. We have used two probes: 74 for the vector (charge) and 1 for the scalar (matter) 
distribution. The radial distributions, x* J (R), are then extracted by fitting the C3 with 

N max 

C,(R,T)^ £ cfiT m to#(R)*r m *c h (2.2) 

where the to,-, c; are those extracted from C2. The most recent calculations are the P+ and DH — 
distributions (see Figs. §-§). Earlier S-wave distribution calculations have been published in 
Refs. §0]. 




3. A model based on the Dirac equation 

A simple model based on the Dirac equation is used to try to describe the lattice data. Since 
the mass of the heavy quark is infinite we have essentially a one-body problem. The potential in 
the Dirac equation has a linearly rising scalar part, &sc^> as well as a vector part by^cR. The OGE 
potential, incorporating the running coupling constant a s {k ), is obtained by the replacement 

r — ► a OGE~ I dkj (kR)a s (k 2 ), a s (k 2 ) = — In" 1 g . (3.1) 

Here Aq£tj = 260 MeV, the dynamical gluon mass m g = 290 MeV (see Ref. [Q] for details) 
and tfQQE is an overall adjustable parameter. The potential has also a scalar term mcoh(h + 1) 
depending on the orbital angular momentum. 
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Figure 3: S-wave ground state (x 11 ) charge distribution. 
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Figure 4: S-wave charge distribution: the off-diagonal correlation between the ground state and the first 
radially excited state. The Dirac model seems to get the node right - while no attempt has been made to fit 
the radial distributions. 
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Figure 5: P-wave charge distributions. 
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Figure 6: P-wave matter distribution. Note that the Dirac model is able to get the sign change right. 
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Figure 7: Comparison of static (infinitely heavy) and smeared ("sum6") heavy quark charge distributions 
for the P+ state. 
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Figure 8: Comparison of static (infinitely heavy) and smeared ("sum6") heavy quark charge distributions 
for the DH — state. 
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Figure 9: Charge sum rule. This should be one (the charge of the q) in the normalisation that is used here. 

The solid lines in Figs. are radial distributions from the Dirac model fit. The fit parameters 
are m = 0.088 GeV, aqgE = °- 81 ' = 1.14 GeV/fm, b vec = 1.12 GeV/fm and co = 0.028. At 
present we only fit energy differences — the ground state energies and the 2S in Fig. || — and sim- 
ply live with the resulting wavefunctions. An attempt could be made to fit the radial distributions 
as well. 

4. Conclusions and acknowledgments 

• The spin-orbit splitting is small and supports the symmetry &vec = ^sc as proposed in 
Ref. [§. 

• There is now an abundance of lattice data to comprehend. The energies and radial distri- 
butions of S, P and DH — states can be qualitatively understood by using a Dirac equation 
model. 

J. K. and A. M. G. wish to thank the UKQCD Collaboration for providing the lattice config- 
urations and the CSC - Scientific Computing Ltd. for providing the computer resources. J. K. and 
A. M. G. acknowledge support by the Academy of Finland (contract 54038), the Magnus Ehrnrooth 
foundation and the EU grant HPRN-CT-2002-003 1 1 Euridice. 

References 

[1] UKQCD Collaboration, A. M. Green, J. Koponen, C. Michael and P. Pennanen, Phys. Rev. D 65, 
014512(2002) 

[2] UKQCD Collaboration, A. M. Green, J. Koponen, C. Michael and P. Pennanen, Eur. Phys. J. C 28, 79 
(2003) 

[3] UKQCD Collaboration, A. M. Green, J. Koponen, C. McNeile, C. Michael and G. Thompson, Phys. 
Rev. D 69, 094505 (2004) 

[4] T. A. Lahde, Nucl. Phys.A714, 183 (2003), hep-ph/02 08110 

[5] P. R. Page, T. Goldman and J. N. Ginocchio, Phys. Rev. Lett. 86, 204 (2001), hep-ph/ 02 94 



T 




ir 



-_ OODF3 

- Q-QDF3sum6 



205/6 



